function [V,phi,u] = update_mesh(V,T,VT,ring1,NodeFlag,phi,u)
% function [V,phi,u] = update_mesh(V,T,VT,ring1,NodeFlag,phi,u)

% update the mesh using the information of phi
% and update the values of phi and u to the new mesh